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Abstract 

If the lepton mixing angle #23 is not maximal, there arises a problem of ambiguity in determining 
#23 due to the existence of two degenerate solutions, one in the first and the other in the second 
octant. We discuss an experimental strategy for resolving the #23 octant degeneracy by combining 
reactor measurement of #13 with accelerator Vu disappearance and v e appearance experiments. 
The robustness of the #23 degeneracy and the difficulty in lifting it only by accelerator experiments 
with conventional Un (and z?«) beam are demonstrated by analytical and numerical treatments. Our 
method offers a way to overcome the difficulty and can resolve the degeneracy between solutions 
sin 2 6*23 = 0.4 and sin 2 9 23 = 0.6 if sin 2 26>i 3 > 0.05 at 95% CL by assuming the T2K phase II 
experiment and a reactor measurement with an exposure of 10 GW-kt-yr. The dependence of the 
resolving power of the octant degeneracy on the systematic errors of reactor experiments is also 
examined. 
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I. INTRODUCTION 



The atmospheric neutrino observation by Super-Kamiokande (SK), which discovered neu- 
trino oscillation, established that the mixing angle #23 is large, which may be even close to 
the maximum |l[. Then, it became clear that the solar mixing angle 612 is also large, though 
not maximal, as indicated by accumulation of the solar and the reactor data 0, 0| • Curiously 
enough, the third angle #13 is known to be small Comprehending the structure of lep- 

ton mixing described by the Maki-Nakagawa-Sakata (MNS) matrix [6| which is composed 
of a nearly maximal, a large, and a small angle should reveal the key to understand the 
underlying physics of lepton flavor mixing. 

Toward the explanation of the coexistence of nearly maximal and small mixing angles 
some symmetries [7| which are motivated by more phenomenological \i <-> r exchange symme- 
tries [8| are discussed. Interestingly, most of them share an attractive feature that #23 = vr/4 
and #i3 = in the symmetry limit, suggesting a possible common origin of the two small 
quantities in the lepton flavor mixing, 13 and -D23 = \ — sin 2 #23 > the parameter which 
measures deviation from the maximal mixing. The quark- lepton complementarity 0, if 
extends to the 2-3 sector, might be relevant for maximal or nearly maximal #23 . If these 
small quantities take non- vanishing values it should give us a further hint for deciding if the 
symmetries are the right explanation for the two small quantities, and if so, for identifying 
the correct symmetry. 

In this paper, we describe an experimental strategy for determining #23 following (lp| . 
where it was proposed that the #23 octant degeneracy can be lifted by combining reactor 
measurement of #13 with accelerator appearance and disappearance measurement of certain 
combinations of #23 and #13. See Ref. |ll| for an earlier suggestion of this possibility. We thor- 
oughly examine this method by taking a concrete (and probably the best thinkable) setting 
for accelerator experiments, phase II of the T2K experiment |l2|, with Hyper-Kamiokande 
as a detector and a 4 MW neutrino beam from the J-PARC facility. For knowledges 
of reactor measurement, we are benefited by the outcome of the world-wide effort We 
carry out a detailed quantitative analysis to obtain the region of #13 and #23 in which the 
6*23 octant degeneracy can be resolved by our method. For a related analysis, see [l5| . 

The #23 octant degeneracy is a part of a larger structure called the parameter degeneracy. 
It may be understood as the intrinsic degeneracy 0] duplicated by the unknown sign of 
Amo i Il7l ] and 6*23 octant ambiguity which lead to the total eight-fold degeneracy fl8| . 
(See (l9| for exposition of analytic structure of the degeneracies.) Then, one has to address 
a problem; Is it possible to resolve only the #23 degeneracy, leaving the other two types 
of degeneracies untouched? We will answer in the positive this question in the context 
of conventional beam experiments. In short, in modest baseline distances for which 
perturbative treatment of the matter effect is valid, the #23 degeneracy is decoupled from 
the other degeneracies and can be resolved independently in the presence of them. 

The current bounds on these small quantities in the lepton flavor mixing are rather mild, 

-0.14 < D 23 = \ -sin 2 #23 < 0.14, (1) 

at 90% CL pj, which is nothing but a translation of the bound sin 2 2# 2 3 > 0.92. At present, 
there is neither indication of deviation from the maximal #23, nor preference of the particular 
octant (apart from very slight preference of the 2nd octant) in the analysis by the SK group 
with their current data set |20J]. On the other hand, the bound on #13 is much stronger if 
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one uses the same variable, 

sin 2 6 n < 0.022 (0.047) (2) 



at 90% (3cr ) CL for 1 degree of freedom, as obtained by the global analysis |2l| with use 
of all the data including the Chooz, the atmospheric, and the K2K |22j one. As discussed 
in (see also one of the major difficulties for accurate measurement of 6*23 with 

accelerator neutrinos is the 6*23 degeneracy. Hence, we expect that our method will help to 
improve the situation. 

In Sec. |nj we explain how the 623 octant degeneracy can be resolved by our method. In 
Sec. IIIII we discuss how robust is the 623 octant degeneracy by indicating the difficulties in 
resolving it only by accelerator measurement. In Sec. II VI we fully explain the statistical 
procedure of our analysis. In Sec. |V| we present the results of our analysis. In Sec. I Vll 
we give concluding remarks. In Appendix [XJ we describe some details of how the referred 
numbers of events are computed in our paper. 



II. THE METHOD AND WHAT IS NEW? 



For completeness of the presentation, we start by reviewing the method for solving 623 
octant degeneracy described in 1(| . At the end of this section, we will try to elucidate 



the difference between this work and the reference [10J. We invite the readers to look at 
Fig. [TJ and first focus on the upper four panels, the case where the input value of 623 
is in the first octant. (The upper four and the lower four panels of Fig. are for input 
values of S23 = 0.458 and 0.542, respectively.) Fig. [T^ describes the constraints imposed 
by each accelerator experiment, u^ (and u^) disappearance and v e (and u e ) appearance 
measurement. Although these contours come from our full analysis whose details will be 
explained in Sees. |TV| and |V| the main features of Fig. [T] can be understood by the vacuum 
oscillation approximation, and essentially it is all that we need. 

The u,, disappearance and u e appearance probabilities in one Am 2 dominance approxi- 
mation |25j, which may be justified by Am^/ Am^ ~ 1/30 -C 1, are given by 

PW - v,) = 1 - sin 2 20 23 sin 2 (^f^) , (3) 
- u e ) = s 2 3 sin 2 20 13 sin 2 (^^) , (4) 

where E denotes the neutrino energy and L is the baseline distance. We use the standard 
notation of the MNS matrix including the symbol Sjj for sin6V,- [2^. Am 2 ^ is defined as 
Amji = m 2 - — mf by using the neutrino masses rrii (i = 1-3). In this approximation P(v a — > 
Up) = P(u a —>■ Up). Accelerator disappearance measurement is expected to determine both 
Amjj and sin 2 2 #23 with high accuracies. It is obvious that, if 623 7^ vr/4, one has two-fold 

solutions of #23 for a given value of sin 2 26*23; s| 3 = ^ 1 ± \f\ — sin 2 26* 2 3 . This is a simple 

illustration of how 6*23 octant degeneracy arises. The almost horizontal two lines in Fig. ^1 
are nothing but these two solutions of S23 for a given value of sin 2 26*23 — 0.993. 

On the other hand, accelerator appearance measurement determines a particular com- 
bination of two angles, S23 sin 2 26*13, as seen in (j3J) for a given value of Am^. The latter 
quantity is expected to be well determined by disappearance measurement. The curved 
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In appearance mode! Am^=2.5 X 10 3 eV J is fixed 
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FIG. 1: The upper (lower) four panels describe the process of how the #23 octant degeneracy can be 
resolved for the case where the true value of sin 2 #23 = 0.458 (0.542), corresponding to sin 2 2#23 = 
0.993. The other input mixing parameters are given as Am^ = 2.5 x 10 -3 eV 2 , sin 2 2#i3 = 0.1 and 
5 = 0, Am 2 ! = 8.0 x 10~ 5 eV 2 , sin 2 612 = 0.31 (the input values of sin 2 26*13 and sin 2 623 are indicated 
by the symbol of star in the plot), (a) The regions enclosed by the solid and the dashed curves 
are allowed regions only by the results of appearance and disappearance accelerator measurement, 
respectively, (b) The regions that remain allowed when results of appearance and disappearance 
measurement are combined, (c) The regions allowed by reactor measurement, (d) The regions 
allowed after combining the results of appearance and disappearance accelerator experiments with 
the reactor measurement. The exposures for accelerator are assumed to be 2 (6) years of neutrino 
(anti-neutrino) running with 4 MW beam power with Hyper-Kamionande whose fiducial volume 
is 0.54 Mt, whereas for the reactor we assume an exposure of 10 GW-kt-yr. The case of optimistic 
systematic error is taken. (See Sees. HVI and Appendix^for details.) 



strip in Fig. [IJi, which takes the approximate form s 2 3 sin 2 26 ) 13 = constant, is the outcome 
of the u e and u e appearance measurement. In Fig. we combine the disappearance and 
appearance measurement and we end up with a pair of degenerate solutions of 613 and 623 ■ 
Now, we combine the information from reactor experiments. In a very good approxima- 
tion what they measure is equal to the vacuum oscillation probability P{v e — > P e ) which is 
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given again in the one Am 2 dominance approximation as 



P(z/ e — * z/ e ) = 1 — sin 2 2# 13 sin 2 



( 



Am 2 .,!, 



) 



(5) 



AE 



Given the value of Am^ by disappearance measurement, the reactor experiments determine 
#13 independent of other mixing parameters. This is illustrated in Fig.^. When we combine 
the accelerator disappearance and appearance experiments as well as reactor measurement, 
one of the two allowed solutions in Fig. [TJo disappears as shown in Fig. [T]i. Thus, the 823 
octant degeneracy can be lifted. 

Though using the same method as proposed in ^(j , our analysis goes beyond that given 
in the reference in a number of ways. We have taken into account errors in accelerator 
disappearance and appearance as well as background so that the event number distributions 
roughly reproduce those obtained by the experimental group |27|, |28(. Our treatment of the 
reactor measurement of #13 is elaborated to include uncorrelated and correlated errors in 
order to treat the so called phase II type high statistics measurement. By using the highest 
sensitivity reactor and accelerator experiments, we aim at revealing the ultimate sensitivity 
for resolution of the 823 octant degeneracy achievable by this method. Differences between 
this paper and jj.lO] exist not only in the method for the analysis but also in some features 
of the results. For the particular set of parameters used in Fig. the resolving power of the 
823 degeneracy is greater for the case with true values of 823 in the first octant than for the 
case in the second octant. It is in disagreement with the naive expectation in [lCj based on 
the difference between 6*23 of a true and a fake solutions. See Sec. |V|for more details. 

III. THE 023 OCTANT DEGENERACY IS HARD TO RESOLVE BY ACCELER- 
ATOR EXPERIMENTS 

We illustrate how difficult is to resolve the octant degeneracy by accelerator experiments 
with conventional neutrino beam. Apparently, this fact has been recognized by people in 
the community, but to our knowledge, a coherent discussion of why it is so has never been 
given. Therefore, we try to fill the gap. As a byproduct, the argument illuminates the 
nature of the 823 degeneracy and may indicate a unique feature of our approach. We note 
that, strictly speaking, our following discussion in this section is valid under circumstances 
that the matter effect can be treated perturbatively, and hence for baseline shorter than 
L ~ 1000 km. (See below.) 

A. The octant degeneracy is robust 

Robustness of the degeneracy is obvious from Fig.[T^ and Fig.[TjD, but we want to elaborate 
the discussion to bring our understanding to a little deeper level. The reasons for the 
robustness are mainly twofold: 

• The difference in energy spectra predicted by degenerate solutions in either the ap- 
pearance or disappearance channels is too small to distinguish between the first- and 
the second-octant solutions. 

• The matter effects in the disappearance channels are not strong enough to lift the 
octant degeneracy, unless one goes to very long baseline. 
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To explain the first point, we introduce the quantity APi 2 (z/q, — > v@), 
AP 12 (z/ Q -> up) = P{y a -> ^; 0*f, 0jf) - P(i/ a -> ^; ft 



2nd fl2nd\ 
23 J P 13 J 



(6) 



the difference between probabilities with parameters of two degenerate solutions, 9\ 3 and 
6\ 3 (i=lst, 2nd), determined at a particular value of energy in this exercise. In Fig. 



plotted are AP 12 (^ AJ 



(left upper panel), AP 12 (z^ 



(right upper panel), and their 



antineutrino counterparts in the lower two panels. As one can see in Fig. |21 AP 12 's are less 
than ~ 0.2% in most of the energy region. In this exercise, the CP phase 5 1 (i=lst, 2nd) is 
kept equal and we did not try to adjust it to make APi 2 smaller. We have also examined 
three other values of delta, 5 = h/2,tt,3tt/2, in addition to the 5 = case in Fig. El and 
reached the same conclusion. 

The smallness of APi 2 (z^ — > u e ) and APi 2 (z/ At — > u^) discourages the possibility of resolu- 
tion of the 23 octant degeneracy by using spectrum informations. This is in sharp contrast 
to the case of intrinsic degeneracy of 613-8, for which the spectrum analysis is proved to be 
powerful in the presently discussed original T2K setting provided that 13 is not too small 



B. Approximate analytic treatment of the # 23 octant degeneracy 



We present an analytic framework to understand better the reasons for the robustness of 
the 623 octant degeneracy. We restrict ourselves to baselines where the earth matter effect 
can be treated as perturbation in the appearance channel. For a consistent treatment we 
start from the expression of the disappearance probability Piv^ — ► v^) in the one Am 2 
dominant vacuum oscillation approximation but with leading order s 2 3 correction 1 



P(v„ 



[sin 2 20 23 + 4s 2 3 s 2 3 (2s 2 23 - l)] sin 



13 

AE 



(?) 



Then, a disappearance measurement determines s 23 to first order in s 2 3 as (s 



>23) 



S23) (1 + s i3) where (S23) is the solution obtained by ignoring s 2 3 . By reading off 
the coefficient of sine squared term in (J7J), one can determine sin 2 20 23 and s 23 is given by 

1 ± \fl — sin 2 26 23 as noted before. We note that the linear dependence of 



'23 J 



1(0) 



( s 23) 011 s i3 is clearly seen in Fig. EH- 

For the appearance channel, we use the v e (u e ) appearance probability with first-order 
matter effect 130] 



v e {v e )\ 
sin 2 26>i 3 s 23 



sin 
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AE 
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=F sm 



Am 31 L 



2P 



2P 



n . . 2 / Am 31 L 

=F 2 sm 5 sm z — 

V AE 



(8) 



1 Notice that the solar oscillation effect which is ignored in (J7J is small because of the suppression by a 
factor of ~ (Am^/Am 2 ^) 2 at around the first oscillation maximum. 
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FIG. 2: Two examples of the #23 degenerate solutions are shown to have almost the same energy 
spectrum in Kamioka with baseline of 295 km from J-PARC. In the upper (lower) panels, we show 
the difference of the probabilities of two degenerate solutions for neutrino (anti-neutrino) channels. 
The solid and the dashed lines are for normal (Am^ > 0) and the inverted (Am^ < 0) mass 
hierarchies, respectively. 



In (JHJ), a(x) = y/2GpN e (x) [31] where Gf is the Fermi constant, N e (x) denotes the electron 
number density at the point x in the earth, J r (= Ci 2 Si 2 c 2 3 Si3C 2 3S 2 3) denotes the reduced 
Jarlskog factor, and the upper and the lower sign ± refer to the neutrino and anti-neutrino 
channels, respectively. 

We make an approximation of ignoring terms of order (Am^i/ Am 31 )J r cos 2# 2 3 • Note 
that keeping only the leading order in this quantity is reasonable because J r < 0.04, 
Am^/Am^ ~ 1/30, and cos2# 23 = ±0.2 for sin 2 26» 23 = 0.96. Then, the two degener- 
ate solutions obey an approximate relationship 

(sin*29 13 sl 3 ) lst =(sin 2 29 13 st 3 ) 2nd , (9) 

or, slfs^f = sl 3 d S2 3 d ignoring higher order terms in si 3 . We can neglect the leading order 
correction in s 2 3 to s 23 in these relations because it gives 0(sf 3 ) terms. 
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C. Appearance channel; v u 



With our machinery well oiled we can understand better the behavior of IS.P^iv^ — > u e ) 
and AP 12 (i^ — > v^) given in Fig. |21 in particular their small values, ~ 10~ 3 . We also note 

and 

AP V 



that even if we switch off the matter effect, there is no visible change in AP 12 (i^ — * v e 

v e ) in Fig. |21 Now, these features can be understood in our analytic framework. 

the difference between 



Noting that J r lst — J 2na! = cos29^ t J^ st in leading order in cos 2# 23 , 
probabilities with the first and the second octant solutions can be given by 



APi 2 (^ -> v e ) 
1,1. co 2U Xi [ ^ E 



. Am^L\ . 2 / Am 2 31 L 

cos <5 sin I — I =F2sm5sm / 1 



2E 



AE 



(10) 



The size of APi 2 is about ~ 10~ 3 for cos26* 23 = 0.2, which is roughly consistent with the 
results given in Fig. El The remarkable feature of (fTUj) is that the leading-order matter 
effect terms drops out completely, because they depend upon 6* 23 and #i 3 only though the 
invariant combination Q. Therefore, we suspect that our statement about robustness of the 
823 octant degeneracy may apply not only to the T2K experiment but also to experiments 
with longer baseline, including the NOz/A project 3^ • 



D. Disappearance channel; —* 

Does the disappearance channel z/ M — > help? Using the approximate relationship © 
and the first order (in s 2 3 ) corrected formula for s 23 mentioned after (J7|), we obtain the 
expression of AP 12 (z/ M — > z/ M ) to first order in matter effect as a sum of the vacuum and the 
matter effect contributions, AP i2 (z/ M — >■ z/ M ) = AP i2 (z/ M — > ^) vac + AP i2 (z/ M -» z^) mat ter- The 
vacuum term 



APi 2 (z/„ -> v. 



fi J vac 



2cos20g' 



- (44)"' + y&gJ?«.s] (^M) sin (^M) (11) 



is small, ~ 10 3 , because of the suppression by either one of s 2 3 or J r , and Am^/Am^. 
The matter term is given by [30( 

APi 2 (z^ -> t' At ) ma tter = 2P(z/ M -> z/ M ) mattor = 2(aL) sin 2 26»i 3 s 2 3 P> 23 P(a;), (12) 

where P 23 = ~ — s 23 , x = A ™g L and the function P is defined as P(x) = | sin 2 (|) — sinx. 
Since the term flips sign under the transformation 623 — > 7r/2 — # 23 , it may be used to 
discriminate if # 23 < vt/4 or 6> 23 > 7r /4 Unfortunately, the size of (fT^ and hence its 
contribution to APi 2 (^ At — > u^) is small, 

sin 2 2# L3 \ / P 23 \ / p \ ( L 



(4P]!U = , xrtH ^ ^ (_£_, (13) 

Note that P is a monotonically increasing function of x in < x < 4, and P(vt) = 4/vt = 
1.273. Thus, the difference (APi 2 ) mat t r is of order 10~ 3 at L = 300 km (less than 1% even 
at L = 1000 km) with #13 at around the Chooz limit 
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It should be noticed that AP i2 flips sign not only under P 2 3 — > — D23 but also under 
Amgj — > -Amjj. Therefore, measurement of the sign of APi 2 determines the combined 
sign, Am^ x P 23 . Hence, the use of this effect to resolve the 823 degeneracy requires prior 
knowledge of the neutrino mass hierarchy, i.e., the sign of Am^. 

Thus, AP 12 is small, typically 0.1% level in the T2K setting, both in the appearance and 
the disappearance channels. This feature explains well the behavior shown in Fig. |21 

E. Approximate decoupling of the 823 octant degeneracy 

Under the approximation of ignoring AP 12 in the appearance and the disappearance 
channels, there is no way to resolve the 623 degeneracy by spectrum analysis. But, on 
the other hand, it implies that the octant degeneracy decouples from the other types of 
degeneracies, Am^-sign and the intrinsic ones. The desirable feature prevails after the 623 
degeneracy is resolved because it is executed by combining the reactor measurement of #13 in 
our method, which is free from any degeneracies. We cannot lift the Am^-sign degeneracy 
by the method we explore in this paper, but we can resolve the 823 octant degeneracy 
independent of the sign of Arafp We will explicitly verify this point in Sec.|Sby performing 
the analysis under assumptions of the right and the wrong mass hierarchies. 

A simple remark on the intrinsic degeneracy of 613-6; The coupling between the 823 and 
the intrinsic degeneracy can be avoided by doin g m easurement at the oscillation maximum, 
or more precisely the "thinnest ellipse" limit [34j. Or, if necessary, it can be resolved 
relatively easily (compared to the Am ', -si gn degeneracy) by doing the spectrum analysis, 
as demonstrated in the T2K setting in [29j . 

IV. ANALYSIS METHOD 

In this section, we summarize the statistical method and the procedure of our analysis. 
We ask the readers to refer to Appenix [X] for any details of how the numbers of events 
are computed. We use the full three-flavor oscillation formulas in our analysis. They are 
obtained by numerically solving the neutrino evolution equation with the constant electron 
number density. 

For concreteness, we consider for accelerator experiment the phase II of the T2K 
project [lH where the beam power is upgraded to 4 MW and the far detector will be 
Hyper-Kamiokande with 0.54 Mt fiducial volume. We consider the 2.5 degree off axis 
beam which has a peak at around 0.65 GeV. We assume the exposures of 2 and 6 years 
of neutrino and anti-neutrino running, respectively [12j. For the reactor experiment, we 
consider the exposure of 10 GW-kt-yr. 

A. — ► u e {p^ — * v e ) appearance mode 

The — > v e {pp —>■ u e ) appearance mode is important to determine precisely the value of 
sin 2 823 sin 2 29 13, provided that A m,, is well determined, which is possible by the disappear- 
ance mode. As discussed in Ref. |34|, if the experiment is done at or close to the oscillation 
maximum, we can determine well the quantity sin 2 823 sin 2 2813, even if we do not know the 
value of the CP phase 6. Since information on the energy dependence of this channel cannot 
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be important in resolving #23 degeneracy (as we saw in the previous section), we consider 
only the total number of events and define x 2 f° r t ne appearance channel as follows, 

I /yobs 1 /yobs /ytheo ytheo \ 2 

2 _ l iv sig "T JV BG iv sig BG I /-. .\ 

Xapp - ^obs + N g + ((Tsigi Vobs )2 + ( aBG iVobs)2 > ^ ) 

where N obs and N theo are the number of events to be observed and the theoretically expected 
one, respectively, for given values of the oscillation parameters. We note that background 
events come mainly from neutral current interactions as well as the events induced by v e (z/ e ) 
which is inevitably contained in the initial flux. We assume optimistic systematic errors, 
<7 S ig = cr BG = 2 % for the T2K II experiment. 



TABLE I: Number of events in the appearance mode with and without oscillation. We assume 
2 (6) years of neutrino (anti-neutrino) running with the T2K phase II set up. For the case with 
oscillation, we set Am^ = 2.5 x 10~ 3 eV 2 , Am^ = 0, sin 2 26*23 = 1> and sin 2 2#i3 = 0.1, with the 
matter density p = 2.3 g/cm 3 and the electron fraction (number of electron per nucleon) being 
Y" e =0.5. Numbers in parentheses correspond to the case where the matter effect is switched off. 



Case 1 (vp — ► v e ) 


Signal events 


BG NC events 


BG beam events 


No Oscillation 





542 


785 


Oscillation 


8683 (8016) 


542 


724 (728) 


Case 2 (p^ — > P e ) 








No Oscillation 





624 


817 


Oscillation 


7340 (7990) 


624 


761 (757) 



In Table I, we show the expected number of events for the T2K phase II for the 2 (6) 
years of exposure for neutrino (anti- neutrino) running with and without oscillation effect. 
For the case with oscillation, we assumed the oscillation parameters Am 2 x = 2.5 x 10~ 3 eV 2 , 
Am^ = 0, sin 2 26*23 — 1> an d sin 2 2#i3 = 0.1. The choice is to compare the number of events 
to the one quoted in Ref. j2?| for T2K I after properly scaling the fiducial volume, exposure 
time and the beam power. We have confirmed that our results agree reasonably well with 
the numbers quoted in [27^ . 

Although we do not use the spectrum informations in our analysis, we present in Fig. |H] 
some examples of the energy distribution of the v e appearance events for completeness. 
For this plot, we set the oscillation parameter Am^ = 2.5 x 10~ 3 eV 2 , sin 2 #23 = 0.4, 
sin 2 2^13 = 0.1, Am 2 ! = 8.0 x 10" 5 eV 2 , sin 2 12 = 0.31, 5 = 0. Throughout our analysis, 
we fix the solar neutrino mixing parameters Am^i and sin 2 6*12 to these values. One can 
observe in the figure that background is quite small in size and has similar shape as the 
signal events. Notice also that the modulation of energy spectrum by neutrino oscillation is 
rather modest. 



B. — > Un (z^t — ► z^t) disappearance mode 

The z/ M — > [p^ — > Vp) disappearance mode is important to determine sin 2 26*23 as well 
as Am| x accurately. In Fig. HI we show the expected event number distribution as a function 
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FIG. 3: Examples of event number distributions are plotted as a function of the reconstructed 
neutrino energy for the appearance mode with and without oscillations. In the latter, the following 
values of the oscillation parameters are used: Am^ = 2.5 x 10 -3 eV 2 , sin 2 #23 = 0.4, sin 2 2#i3 = 0.1, 
8.0 x 10~ 5 eV 2 sin 2 12 = 0.31, 5 = 0. 



of the reconstructed neutrino energy in the absence (left panel) and in the presence (right 
panel) of oscillation with the mixing parameters, sin 2 #23 = 0.4, Am 2 3 = 2.5 x 10~ 3 eV 2 , 
sin 2 26* 13 = 0.1, and 6 = 0. Unlike the case of the appearance channel, the energy distribution 
is significantly modified by the oscillation effect. 

For the disappearance mode, we consider 36 bins with 50 MeV width from 0.2 GeV to 2.0 
GeV in terms of the reconstructed neutrino energy. The \ 2 function is denned as follows, 



2 = . ^ l^f bs + ^ S G - (1 + a sig )N^° - (1 + a BG )flg%] 2 (a si X (a BG \ 2 

^ ~ oSSo V + ^°BG V J + Wg J ' 

(15) 

where and A^ e ° are, the number of signal events to be observed and the theoreti- 

cally expected one, respectively for the z-th bin, and N°]? i and N^r° are the corresponding 
background event numbers. We assume the optimistic values for the systematic errors, 
c s ig = o"bg = 2% for the T2K phase II. 

In FigJSJ we show the expected sensitivity for sin 2 2^23 assuming the pure 2 flavor oscil- 
lation (6 , 13 = 0) as a function of the true value of Am^, which compares reasonably well 
with that of j28|. We have checked that this simple x 2 can reproduce the expected T2K 
sensitivity based on more realistic calculations obtained by Monte Carlo simulation given in 
Ref. @. 
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FIG. 4: Examples of event number distributions are plotted as a function of the reconstructed 
neutrino energy for the disappearance mode without oscillation (left panel) and with oscillation 

2.5 x 10~ 3 eV 2 , sin 2 23 = 0.4, sin 2 20 13 = 0.1, 



Anij! = 8.0 x 10 5 eV 2 , sin 2 #12 = 0.31 and 5 = 0. The histogram by the solid lines indicate the 



(right panel) with the mixing parameters Am^ 

sin 2 #12 

sum of the signal and background events whereas the ones by the dashed line indicate only the 
background events. 
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C. v e — > v e disappearance mode 



For the reactor experiment, we use the same x 2 function used in our previous works 
36[, which is defined as, 



Xreac = m j n Y\ 
as ' * 

a=f,n 



f (N£~>- + * + + a)N$*) 2 a 2 { ,2 
^ 1 N$* + ^ b (^ 8 ) 2 <*u I cr 



a 2 , , 
+ -5-,(16) 

°DB 



where iV^ eo represents the theoretical number of events at a near (a = n) or a far (a = /) 
detector within the z-th bin whose width is 0.425 MeV. We set the distance to the near and 
the far detectors to be 300 m and 1500 m, respectively. N°^ s are number of events to be 
observed. We consider four types of systematic error: o"db, 0"Db, 0"dB and crab- The subscript 
D (d) represents the fact that the error is correlated (uncorrelated) between detectors. The 
subscript B (b) represents the fact that the error is correlated (uncorrelated) among bins. 
The a's are parameters to be varied freely in order to take into account these systematic 
errors. In this work, we consider the high precision reactor experiments whose sensitivity 
can go beyond the ones currently expected by experiments such as Double-Chooz [3?} and 
KASKA j38|. Namely, we assume sensitivities below sin 2 26> 13 = 0.01, the one expected to 
be achievable b y p hase-H type experiments such as the Braidwood [3£|, the Daya Bay |4pj . 
and the Angra |4l| projects. 

Computation of the event number is done in the same way as in Ref . j36| . We ignore the 
possible contribution from geo-neutrinos. It is demonstrated recently by KamLAND that 
its flux is consistent with the one expected by geo-chemical earth models j4^|. In this case, 
the effect of geo-neutrinos is negligibly small at the baseline of ~1 km, as one can easily 
guess by extrapolation of the situation in a reactor # 12 experiment with baseline of ~60 km 

To characterize reactor measurement, we use GW-kt-yr, the "total exposure" unit, which 
is defined as the product of the net values of the reactor thermal power (in GW), detector 
fiducial volume (in kton) and running time (in year). We consider the reactor measurement 
for 10 GW-kt-yr. The total number of events at the far detector is 1.63 x 10 6 . We consider 
two different sets of systematic errors: a relatively conservative choice and an optimistic 
one. For the conservative choice, we adopt similar values of the systematic errors used in 
Ref. 35l 36|, o"db = CDb = 2.0 % and a dB = a dh = 0.5 %, and for the optimistic one, we set 



odb = °Db = 1.0 % and <j dB = 0.2% and cr db = 0.2 %. The latter extremely small errors may 
be difficult to reach, but they are used to estimate the upper limit of resolving power of the 
623 degeneracy by the present method. The sensitivity limit of #13 (assuming no depletion) 
at Am§! = 2.5 x 10" 3 eV 2 is sin 2 26 n = 1.14 x 10~ 2 (2.63 x 10~ 2 ) at la (3a) CL for the 
relatively conservative errors, and sin 2 2^13 = 5.39 x 10~ 3 (1.25 x 10~ 2 ) at la (3a) CL for 
the optimistic ones. 



D. Combined analysis 

For the combined analysis, we simply sum all the \ 2 functions defined in Eqs. (J 14)1 . (|15|) 
and ([TBI. 

2 — 2,2,2 (~\7\ 
X Xapp ' Xdis "1 Xreac " \ / 

The allowed region in the sin 2 2#i 3 — sin 2 623 plane is determined by the usual condition, 
Ax 2 = X 2 ~ Xmin < 2.3, 6.18 and 11.83 for 1, 2 and 3 a CL for two degrees of freedom. We 
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will establish the parameter regions where we can resolve the 6*23 octant degeneracy for 1 
degree of freedom by imposing the condition |Xmin($23 U °) — Xmin($23 SC )l > 2.71, 4 and 6.63 
for 90, 95 and 99% CL, respectively, where 6*23 Ue and 6*2| lse are > respectively, the true and the 
false value of 6*23. 

V. ANALYSIS RESULTS 

In this section, we show our results based on our analysis with the combined \ 2 °f all 
channels. To remind the readers, our analysis is based on the T2K II experiment of 2 (6) 
years running of neutrino (anti-neutrino) modes with 4MW beam power with the Hyper- 
Kamiokande detector whose fiducial volume is 0.54 Mt [lil ]. For the reactor experiment, the 
exposure of 10 GW-kt-yr is assumed. 

We assume throughout this section that the mass hierarchy (determined by nature) is 
normal type (Am^ > 0) unless otherwise stated. Even if we repeat the same procedure with 
the true mass hierarchy of inverted type (Am| x < 0), the region in which the degeneracy is 
solved is remarkably similar. Therefore, we decided to concentrate on the normal hierarchy 
case. Of course, we examine the stability of our results by assuming the wrong hierarchy in 
the analysis. 

In Fig. El we show the allowed regions in the sin 2 26* 13 — sin 2 6*23 plane at 1 a (light blue 
curve), 2 a (black curve), and 3 a (red curve) CL for 2 degree of freedom, for 8 different sets 
of input parameters; (sin 2 26* 13 , sin 2 6* 23 ) = (0.08, 0.4), (0.08, 0.6), (0.08, 0.45), (0.08, 0.55), 
(0.04, 0.4), (0.04, 0.6) and (0.04, 0.45), (0.04, 0.55), which are indicated by the symbol of 
a star. They are determined by the combined x 2 analysis using all the channels, — > v e 
{pp — > v e ) appearance mode, by — > (z/ M — > v^) disappearance mode and by v e — > v e 
disappearance mode. The relatively conservative values of the systematic errors are taken 
for the reactor measurement; o"db = ^Db = 2.0 % and a^B = o"db = 0.5 %. In our \ 2 analysis, 
for given values of input parameters, we vary not only 6*23 and 6* 13 but also S and Am^, as 
these parameters should be determined by the fit. Note, however, that the range of Am^ 
and #23 are restricted to the ones constrained by atmospheric neutrino experiments p|. 

For a larger value of sin 2 26* 13 = 0.08, the octant degeneracy is resolved for sin 2 26* 2 3 = 0.96 
(see Fig.^l and bl), whereas for sin 2 26*23 = 0.99, the degeneracy can be resolved only at 1 a 
CL but not at 2 a CL or higher (see Fig.^fcl and dl). For a smaller value of sin 2 26*13 = 0.04, 
for sin 2 2^23 = 0.96, the octant degeneracy is resolved only at 1 a CL (see Fig. E^2 and b2) 
and for sin 2 26*23 — 0.99, the degeneracy can not be resolved even at 1 a CL (see Fig. Efc2 
and d2). 

In Fig. we show the same quantities but for the small systematic errors for reactor 
experiments, a DB = o"Db = 1 %, cr db = a dB = 0.2 %, the very optimistic ones. In this 
case, the octant degeneracy is completely resolved for a large #13, sin 2 2# 13 = 0.08 both 
for sin 2 26*23 = 0.96 and 0.99 (see Fig. [7^,1, bl, cl, and dl). For a small value of #13, 
sin 2 26*13 — 0.04, however, we have a mixed result; For sin 2 26*23 — 0.96, the 6*23 degeneracy 
is resolved both for sin 2 26*23 — 0.96 and 0.99 apart from a tiny 3a CL region in Fig. [7^2. 
For sin 2 26*23 = 0.99, the degeneracy is not lifted, except in Fig. [TJi2 where the clone solution 
disappears at la CL. 

In Fig. IS] and Fig. |UJ we show the region bounded by the solid and the dashed curves in 
which we can resolve the 6*23 octant degeneracy in the sin 2 6*i 3 -sin 2 6*23 plane at 90% (thin 
green curves) and 99% (thick red curves) CL for 1 degree of freedom for the conservative 
and the optimistic values of the reactor systematic errors, respectively. Namely, the region 
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FIG. 6: Some examples of the allowed regions in the sin 2 #13 — sin 2 623 plane determined by 
the combined analysis for various different input parameters. For all plots, the input values of 
Anijj = 2.5 x 10 _3 eV 2 and the CP phase 5 = 0. The input value of #13 for the upper (lower) four 
panels is sin 2 6*23 = 0.08 (0.04). The input values of 6*23 are (a) sin 2 6*23 = 0.4, (b) sin 2 6*23 = 0.6, (c) 
sin 2 623 = 0.45 and (d) sin 2 623 = 0.55. The conservative values of the systematic errors are taken 
for reactor measurement; <tdb = CDb = 2.0 % and o^B = <?db = 0.5 %. 



inside the bands correspond to the parameters which satisfies |Xmm(^23 UC ) ~~ Xmin($23 SC )l > 
2.71 and 6.63 where 9% 3 ue and 6*2| lse are the true and the false values of 6*23, respectively. The 
solid (dashed) curves in Fig. |E]and Fig. are for cases that we fit the data set generated with 
Am| x > under the hypothesis of right hierarchy, Am| x > (wrong hierarchy, Am^ < 0). 
If we are ignorant about the mass hierarchy the case of worse sensitivity (wrong hierarchy) 
must be considered as the sensitivity region. The fact that the solid and the dashed curves 
come close to each others proves that the resolution of the 623 octant degeneracy can be car- 
ried out independent of the lack of knowledge on the neutrino mass hierarchy, the decoupling 
of the two degeneracies as argued in Sec. IIIIE1 
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FIG. 7: Same as Fig.Elbut for smaller systematic errors for the reactor experiment, <tdb = <7Db = 1 
%, a dh = 0.2 % and a dB = 0.2 %. 



The figures indicate that our method of combining reactor measurement of #13 with 
the accelerator disappearance and appearance experiments allows to resolve the 823 octant 
degeneracy to a reasonable level. This is highly nontrivial because it is quite a robust degen- 
eracy which is hard to lift by using only relatively short-baseline (~ 1000 km) accelerator 
measurement, as we saw in Sec. IIHI We can observe by comparison between Fig. |H]and Fig. 
El that improvement of the systematic errors in reactor measurement is the key to the better 
resolving power of the octant degeneracy. It should also be mentioned that the results in 
Fig. 01 are obtained with the extremely small systematic errors. Hence, the results may be 
interpreted as the limit achievable by the present method. 

So far we have assumed that the input mass hierarchy is normal while examining the true 
and the false mass hierarchies in analyzing data. In order to verify that resolving power of 623 
degeneracy is not affected by the unknown sign of Am^, we tried the converse. Namely, we 
tried to fit the data set generated with the inverted hierarchy Am^ < under the hypothesis 
of the right and the wrong hierarchies. We do not show the results because they are very 
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90% (thin green) and 99% (thick red) CL. The solid (dashed) curve is for the case of taking the nor- 
mal (inverted) hierarchy to perform the fit, assuming the normal hierarchy as input. Conservative 
systematic errors, as indicated in the figure, are considered here. 



similar to those for the input normal mass hierarchy. Of course, this is expected because 
the experimental setting we consider in this work is not sensitive to the sign of Am^, which 
in turn implies that our results must be insensitive to the mass hierarchy confusion. This 
point was discussed in depth in Sec. Illll 

We note that globally our results in Fig. |H] and Fig. El are similar to what we can find in 
Fig. 7 of Ref. jldj | . which was obtained by a much more simplified analysis. However, there 
are some notable differences. In Fig. 7 of Ref. [lfjj . it is always true that the 823 degeneracy 
is easier to lift in the second octant, 623 > 7r/4. But, we observe that it is true only for 
relatively small 9 13 , sin 2 2^13 < 0.06. For relatively large 6*13, sin 2 2#i 3 > 0.06, the octant 
degeneracy is easier to be resolved for the first octant case 623 < it /A than for the second 
octant. 

It appears that at large 8 13 , the two effects which were not taken into account in the 
treatment in ^(J come into play, as indicated in Fig. d The width of the appearance strip, 
which comes from ^-dependent terms in the oscillation probability (see Eq. (JHJ) in Sec. IIII B|) . 
is narrower at smaller #i| ke , which is the case of true 623 in the first octant, making rejection 
of the fake solution easier. Also, the two disappearance "lines" come closer in the case of 
true 623 in the second octant, which produces the similar effect. As a result of the two 
effects which simultaneously act toward the same direction, the degeneracy turns out to 
be easier to resolve for the case of true 623 in the first octant. We, however, want to note 
that the feature of first-second octants asymmetry, in particular in Fig. |S1 depends upon the 
systematic errors and its full understanding may require much more subtle discussions. 
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VI. CONCLUDING REMARKS 

In this paper, we have explored the possibility of resolving the 823 octant degeneracy 
by combining reactor measurement of #13 with the possible hig hest accuracy accelerator 
disappearance and appearance measurement, as proposed in It utilizes the nature of 
the reactor experiment as a pure measurement of #13 to resolve the degeneracy. It is nice 
to see that the reactor measurement can contribute to explore the two small quantities in 
lepton flavor mixing, # 13 and a deviation of 823 from the maximal value. The results of our 
quantitative analysis indicates reasonably high performance of the method. As shown in the 
summary figures of the resolving capabilities, Fig. |S] and Fig. the degeneracy is resolved 
in region of relatively large #13 and a sizable deviation of 623 from the maximal. 

Prior to the quantitative analysis based on our method, we have discussed the robustness 
of the 823 octant degeneracy and illuminated the difficulty in resolving it only by accelerator 
experiments. In particular, we have demonstrated by showing a sample figure, Fig. that 
neither the spectral information nor the matter effect enables us to resolve the degeneracy, 
This feature was also understood on the basis of analytic treatment using the approximate 
formulas for appearance and disappearance probabilities valid to first order in matter effect. 
Considering the robustness of the degeneracy, the opportunity of resolution offered by our 
reactor-accelerator method is highly nontrivial. 

Finally, several remarks are in order: 

(1) As indicated in Fig. ^ the sensitivity of our method for resolving the 823 degeneracy 
is limited mainly by the accuracy of #13 determination by reactor experiments. The region 
without resolving power which remains in Fig. |U] must be taken as the intrinsic limitation 
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of our method for resolving the #23 degeneracy in its current form, because we have already 
assumed a rather optimistic values of systematic errors in the reactor measurement, in 
addition to the extreme accuracies of accelerator experiments. 

(2) In this paper, we have considered the setting of T2K experiment originally described 
in ^ would be interesting to examine how the sensitivity of resolution of #23 degen- 
eracy changes if we adopt the Kamioka-Korea identical two-detector setting by which the 
degeneracy related to the neutrino mass hierarchy and the intrinsic one can be resolved |29| . 

(3) Our method for resolving 823 degeneracy is by no means unique. The other possibilities 
include: accelerator measurement with silver channel [y e — > v T ) |43J] which has a different 
623 dependence, and detection of the solar oscillation term by either atmospheric neutrino 
observation, 3> HI 3] or very long baseline accelerator experiments j4?|. Combination of 



the accelerator and the atmospheric neutrino experiments can also be pursuit |48j . The ulti- 
mate possibility would, of course, be the "everything at once" approach which combines 



measurement by superbeam and neutrino factories with golden and silver channels. 



APPENDIX A: CALCULATION OF NUMBER OF EVENTS 

In this appendix, we provide some detailed informations about how the number of events 
are computed. For clarity of notation we denote the neutrino energy (E in the text) as E v 
in this appendix. 

(i) Apparence channel — > v e (or z/ M — > u e ) 

As we mentioned in sec. IV, since the information on the energy spectrum of the appear- 
ance channel is not important in resolving 823 degeneracy, for this channel, we do not try to 
make binning but use the expected total number of events, which is computed as, 

N sig = n N T dE v (j) v ^E v )P{p^P e -E v )aZ t {E u )e v& {E v ), (Al) 



where E v is the (true) neutrino energy, is the number of target nucleons in the detector, 
T is the running (exposure) time, (j) Ufi (E u ) is the flux spectrum of the 2.5 degree off axis 
beam a^ t (E u ) is the total cross section and e Ue (E u ) is the detection efficiency, which is given 
as a function of the neutrino energy. The detection efficiency we used takes into account all 
the cut imposed to reduce background, which enables us to obtain the signal distribution in 
the range between 0.35 and 0.85 GeV in terms of reconstructed neutrino energy |27| . 
We consider two kinds of background events, 

i\rBG _ atBG(NC) , r BG(beam) 
iV total — iV total + iV total J 

where N^^ c ^ comes from neutral current interactions and iV t ^j beam ^ implies the event 
induced by the v e neutrino existed in the original beam. As for the signal, we used the 
efficiency functions which correspond to the cuts found in Ref. jl?} in order to compute the 
background. 



(ii) Disapparence channel — ► (or — > v^) 

For the disappearance mode, it is essential to make binning in order to take into account 
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the information on the energy spectrum. The expected number of signal events for the i-th 
bin is computed as follows 

e; 1 ; max J?max 

iV si g,i = n N T / dE mc / dE v <f> Vlt {E V )P(^ -> ^; E v )a v »{E v )e™(E rec )R{E v , E mc ), 

(A3) 

where E rec implies reconstructed neutrino energy, a^ C Q E (E Iec ) implies charged current quasi 
elastic (CCQE) reaction cross section, e™(E rec ) is the detection efficiency as a function of the 
reconstructed neutrino energy, which corresponds to the one used in Ref. [281] . R(E V , E Tec ) 
is the Gaussian-like resolution function, 



R(E V , E iec ) = — = — exp 



1 / E v — E rcc 

2 V °e 



(A4) 



where we set cr^ = 80 MeV. 

For this mode, we consider two kinds of background, 

N BG = jv BG(NC) + jv BG ( cc " N Q E ) (A5) 

where j\T BG '- NC ' ) is the events coming from neutral current reaction and j\T BG ^ CC N ® E ' ) is 
the one coming from the charged current non-quasi elastic (CC-NQE) reactions. We note 
that jV BG( ' NC ' ) do not depend on oscillation parameters, whereas _/V BG( ' CG NOjE ^ depend on 
oscillation parameters in a non-trivial way. For the background coming from CC-NQE 
reactions, E iec = E tTue is not a good approximation because events induced by higher en- 
ergy neutrinos mimic the signal induced by lower energy ones. However, we observe that 
it is a good approximation to take that E mc — E tIue = 300 MeV or A ? " BG( ' CC ~ NC ^ E ^(£' rec ) ~ 
Af G(CC ~ NQE) (300 MeV + E tiue ). This allows us to compute V BG(CC " NQE) (E rcc ) in the pres- 
ence of oscillation provided that we know the distribution jy BG ( cc_N Q E ) as a f unc tion of 
reconstructed neutrino energy in the absence of oscillation for the 2.5 degree beam. 

(hi) Disappearance channel v e — > v e 



We compute the expected number of v e events in the i-th energy bin, N^ heo = N* ea ' 



where NJ eac 



N v ™ c (sm z 6 12 ,Am z 21 ) = N P T 

x J dE u e<j>{E v )P{p e ^p e ;L,E u )a{E v ) JdEe det R{E,E'), (A6) 

where N p is the number of target protons in the detector fiducial volume, T is the exposure 
time, and 4>{E U ) is the neutrino flux spectrum from the nuclear power plant (NPP) expected 
at its maximal thermal power operation, e denotes the averaged operation efficiency of the 
NPP for a given exposure period and it is taken to be 100% here under the understanding 
that the unit we use GW t h-kt-yr refers the actual thermal power generated, not the maximal 
value. P{p e v e , L, E v ) is the familiar antineutrino survival probability in vacuum, which 
is given by Eq.flHJ) and it explicitly depends on Am^ and sin 2 9 12 . (t(E u ) is the u e absorption 
cross-section on proton, edet = 0.898 is the detector efficiency and R(E, E') is the energy 
resolution function, which is assumed to have a Gaussian form with E (= -E pr ompt) the 
observed prompt energy (total e + energy) and E' — E v - 0.8 MeV the true one. 
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